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HIGH TEMPERATURE EQUATIOH OF STATE FOR AlUMIHUM 1 

By 

R. J. Naumann 


SUMMARY 

An equation of state is developed that is capable 
of describing metallic elements in both solid and 
liquid phase from ambient temperature and pressures 
to states of extreme temperature and pressure. A 
novel technique for treating the atomic vibrational 
contributions is introduced by assuming each atom 
vibrates independently in a tan 2 potential. The 
Schroedinger equation can be solved exactly for this 
potential function, and the quantum mechanical parti- 
tion function is computed directly. In this manner 
a continuous, thermodynamically consistent descrip- 
tion of metals is obtained which reduces to an 
Einstein solid at moderate temperatures and to an 
ideal gas in the limits of high temperatures and/or 
large volumes. 

Detailed numerical computations are carried out 
for aluminum. All the constants required to 
configure the model are obtainable from elementary 
thermodynamic data at ambient conditions. Excellent 
agreement is obtained with experimental shock 
compression data on solid and porous samples at 
pressures to 5 Mb and temperatures to 25 000°K. 

INTRODUCTION 


Substantial progress has been made in recent 
years in understanding the behavior of metals at 
extreme temperatures and pressures. Shock com- 
pression techniques have achieved pressure measure- 
ments to approximately 10 Mb and temperatures of 
tens of thousands degrees [ 1] . At the lower pres- 
sures, the Debye or Einstein model for energy and 
the Mie-Gruneisen equation for pressure is adequate 
and zero-degree isotherms were extracted from shock 
measurements by using these relations to subtract out 


the thermal contributions. It was found that a Morse 
potential predicted the zero-degree isotherms to a 
fair degree of accuracy [2], which confirmed an 
earlier suggestion by Slater [3] . 

As experimental pressures and temperatures 
increased, various inadequacies in the theory became 
evident. Electronic contributions were introduced. 

For compressible metals, anharmonic terms became 
significant. These were treated by Al’tshuler et al. 
using the free volume theory of Lennard -Jones and 
Devonshire [4] . Pa stine [5] use'd the perturbation 
method of Liebfried and Ludwig [6] to correct for 
anharmonic effects. 

For the high -temperature states reached by 
shocking porous samples, it was found that the 
behavior became more ideal gas-like. Kormer [7] . 
proposes a set of empirical interpolation equations 
to transform solid -like behavior of energy, pressure, 
and heat capacity to the ideal gas relations. Urlin [8] 
proposed an empirical free energy function to 
account for phase transitions. All of these empirical 
relations require adjustable constants which are 
evaluated from experimental high-pressure data. 
Furthermore, there is no guarantee of thermody- 
namic consistency between them 2 . 

An equation of state for metals is developed in 
this work thgt is capable of describing the liquid-dense 
vapor phase as well as the solid phase from ambient 
conditions to temperatures and pressures exceeding 
the present experimental range. Rather than use in- 
terpolation equations to transform the behavior of a 
solid to that of a gas at high temperatures or at low 
densities, the approach will be to start with empirical' 
interatomic potential functions and develop the entire 
equation of state using quantum statistical mechanics. 
The only empirical constants required are obtainable 
from elementary thermodynamic quantities such as 


1. This report was based on work submitted in partial fulfillment of the requirements for the degree of 
Doctor of Philosophy in the Department of Physics in the Graduate School of the University of Alabama. 

2. Kormer’ s pressure, energy, and heat capacity were chosen in such a manner to be consistent among 
themselves. 



heat of vaporization, compressibility, entropy, etc. 

In this manner the behavior of a metal throughout the 
entire high-pressure, high-temperature regime can be 
predicted from its properties at ambient conditions. 
Since all the thermodynamic properties are derived 
from the partial function, thermodynamic consistency 
is guaranteed. This allows a complete thermodynamic 
description , so that phase transitions are described 
naturally from Gibbs energy considerations. 


relation 

p= . (9F k /8V) T 


and is 
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FORMULATION 


Vibrational Terms 


It will be assumed that the Helmholtz energy can 
be expressed as the sum of three independent 
contributions: the lattice term, F^, the ionic 

vibration term, F , and the free electronic term, F . 

v e 

Since all thermodynamic functions are derivable 
from the Helmholtz energy, it immediately follows 
that they all can be expressed as sums of these three 
contributions. The fact that there are interactions 
between the free electrons and the lattice will be 
accounted for by introducing an effective electron 
mass. 


At moderate temperatures, the ionic vibrational 
contributions may be computed from the assumption 
that the ions behave as harmonic oscillators. This 
assumption leads to the Einstein model, in which it is 
assumed that all atoms vibrate independently at the 
same frequency, or to the Debye model which considers 
a distribution of normal mode frequencies. At very 
high temperatures, the vibrational amplitudes are 
such that nonlinear restoring forces must be con- 
sidered. This is accounted for by assuming each atom 
vibrates independently 3 in a potential well given by 


2a 2 mu; 2 

*(x) — 5T- 



(3) 


Lattice Terms 


Since the Morse potential has been found to be a 
reasonable empirical representation for the lattice 
energy of a metallic solid, the Helmholtz energy 
for the zero-degree lattice is written 



where L is the zero-degree heat of vaporization, 

° l/3 

or total binding energy, | = (V/V q ) , and b is an 

empirical constant determined from the compres- 
sibility. From the relation F - E - TS, the internal 
energy, E^, of the zero-degree lattice is identical to 

F^. The zero -degree isotherm is found from the 


When the displacement x is small compared to the 
atomic spacing a, the potential reduces to 

ma, E 

$(x) = — x 2 , (4) 

where t 0 is the Einstein frequency. For a harmonic 

E 

oscillator, w 2 = K/m, therefore, equation (4) has the 
form l/2 Kx 2 where K is the spring constant 
corresponding to the Einstein frequency. 

At larger displacements, x«±a, the potential 
approaches infinity which describes a hard sphere 
collision between point masses. This roughly cor- 
responds to nuclear collisions between neighboring 
atoms whose position expectation values are ±a. 

The choice of the tan 2 functional representation 
was made because it behaves as desired in the limits 


3. The assumption of independent vibration is justified because even at moderate temperatures the Einstein 
and Debye models give almost identical results. 


2 


J 



and also because it allows a closed-form -eigenvalue 
solution to the Schroedinger equation. The energy 
eigenvalues are 4 


where y is the Gruneisen ratio, -(8 In co/d In V) . 
Formulation of y for solids has received extensive 
treatment in the literature. 


e n - (n*+n) « o +n(h*c^+«* o )*/* , 

(5) 

n = 0, 1,2, . . . 

where is the degeneracy energy tt 2 # 2 / (8 m a 2 ) . 

The vibrational partition function Z for N 
atoms, each with three degrees of freedom, is 

oo 

In Z = 3N In Z e~V kT (6) 

n=0 

where k is the Boltzmann constant. 


The brackets ^ ^ denote the ensemble average; eg. , 

oo 

z 


-e /kT 
x e n 


<*> = 


n=0 ' 


00 — € /kT 

Z e " 

n=0 


(13) 


At solid densities, JKco » and the Boltzmann 

factor reduces to exp [-(n 2 e Q + nK^)/kT] . Unless 

kT » Xw, the Boltzmann factor will become negligibly 
small at low enough values of n to prevent the n 2 ^ 

term from contributing significantly. In this case, 
the Helmholtz energy reduces to 


The various thermodynamic functions are 
obtained from the partition function in the usual 
manner. 


— e /\c r V 

F = -kT In Z = -3NkT In ^ e n , (7) 

1 n=0 

E v = kTZ HTF = 3N < e n) > < 8 > 

<*-(£),-■& [W> -«/}■ 

(9) 

and 

P v = - kT ^ lnZ = 3N < P n) ’ (10) 

where 


F = - 3NkT 2 In T e 
v n=0 

= 3NkT In (i-e"^ w/kT ) , 


(14) 


which is the well-known result from the Einstein 
model . 


The Liquid Potential Function 


The zero-degree potential energy in the liquid 
phase is represented by 


Lm e 


BV 

2b( 1-1) 2lo 


4 


. {S{i 


(15) 


Differentiating equation ( 5) , the P contribution 
becomes 


P = 
n 



(X 2 w 2 + e 2 ) 1/2 
o 


’( 12 ) 


where L 1} B, E^, and are empirical constants 

peculiar to the liquid phase, and L q and b are the 

same as in the solid potential, equation ( 1) . This 
particular form was chosen for the following reasons. 


4. See problem 12 in D. ter Haar, Pr oblems in Quantum Mechanics, Academic Press, New York, 1960. 

Actually, there should be a l/2 added to n in both terms corresponding to the ground state. However, 

the ground state energy is included in the zero-degree lattice terms; therefore, it is suppressed in the 

vibrational terms. 
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For ^ the repulsive term is the same as in the 
Morse potential, but the attractive term is the form 
resulting from van der Waals forces. This potential 
function cannot be extended to small £ because 
eventually the attractive term will override the 
repulsive term. To avoid this, the form given for 
£ < £ t was chosen. This is identical to the Morse 
curve for the solid except that and £* have 

been introduced to account for the heat of fusion and 
different compressibility of the liquid phase. 


The pressure is again given by - ( QF^/BV) ^ 
and is 


p k = 


BV 


2bL^ Zb(l-S) o 
3£ Z V 


2b L 


0 


3^rv n 




(16) 

; 


Of the four additional constants introduced in this 
potential, two are required to match the F^ and 

at 1 = ^. The remaining two are chosen to produce 
the observed energy and pressure at the melting 
point 5 . Figure i compares the liquid potential with 
the Morse potential for aluminum. 


definite position in the lattice. In a liquid, particles 
can exchange positions and thereby lose their 
distinguishability. The partition function must be 
adjusted accordingly for "proper Boltzmann counting." 
Since each atom occupies a cell with volume (2a) [3], 
a volume containing N atoms spaced at average ^ 
distance, a, can have N/8 cells. There are (N/8) 
ways of arranging N particles among N/8 cells, of 
which N! are redundant because of particle indistin- 
guishability. The partition function becomes 


Z = 



3N £ 
n=0 


-e /kT 
n 


(17) 


The Helmholtz energy becomes, using Stirling’s 
approximation, 


F = NkT (In 8 
v 


- i) -3NkT In £ 
n=0 


-e /kT 
n 


(18) 


The energy eigenvalues, e , are given by 

equation (5) ; the internal energy and heat capacity 
are given by equations (8) and (9), respectively, and 
the pressure is given by equations ( 10) and ( 12) , 
except that the Einstein frequency and the Gruneisen 
ratio will have different values for the liquid. 


The Liquid Vibrational Component 


The Lennard-Jones Devonshire model of liquid 
and dense vapors treats each atom as though it were 
moving in a cage of its nearest neighbors which are 
considered fixed . Since the number of nearest 
neighbors is usually 12, a high degree of symmetry 
exists, and a spherically symmetric potential maybe 
assumed. Following these concepts, the same form 
of potential function as assumed for the solid, 
equation (3) , is used to describe the behavior of an 
atom in such a cage. 

The partition function will be somewhat different 
from equation (6), however, because of particle 
indistinguishability. In the solid phase, particles 
are, in a sense, distinguishable because of their 


It is interesting to examine the behavior of the 
liquid in the limit of large volumes where tfcu-K). 

In this limit, the atom moves freely within its "cage." 
Encounters with the wall correspond to elastic point 
collisions characteristic to an ideal gas. Since 
kT » e^, the summations maybe replaced with 

integrations and 

CO 

Z e 

n=0 

The Helmholtz energy becomes 



(n z + 2n)€ 

o 

kT 


7rkT 

4e 


i/2 


(19) 


5. The parameter, £ t , is the value of £ corresponding to minimum potential in the liquid phase. It will 
generally be more than unity because the distribution of the atoms in liquid metals is random instead of the 
more efficient close packing usually associated with crystal lattices. This choice of functional dependence 
was made to assure P (£*) = 0 and to make the F and P approach the same values as the solid phase 
at small £. 
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Figure i. Comparison of solid and liquid potentials. (The liquid potential approaches the solid at high 

compressions, differs by approximately the heat of fusion at V = V , and exhibits a van der Waals 

o 

behavior at large volumes. ) 


Similarly, the other equations reduce to the ideal gas 
2 NkT 

relations, E NkT, P = — — , and 

v 3 v V 

3 

C w ” ~ 2 ~ Nk * Usin S the result of equation (20) , the 
relations S = (E-F)/T, and the fact that E k = F , 
the entropy in this limit becomes 



which is the well-known Sakur -Tetrode equation for 
the entropy of an ideal gas. 

At liquid densities, » € and the treatment 
o 

of the vibrational component of liquids at moderate 
temperature reduces to the Einstein model. However, 
the Einstein frequency and the Gruneisen ratio will be 
different in the case of a liquid. 


5 



A crude estimate of the vibrational frequency can 
be made by considering an atom at the origin with 
an atom at ±a. The change in potential resulting 
from a displacement x is 


$(x) = $(a+x) + $(a-x) - 2$ (a) (22) 

Using only the repulsive term in the Morse Potential, 
2b ( i-r/a ) 

$(r) = $ e ° , 

o 


$(x) = e 2b(1 ^ 

o 

For small displacements, 




. 2b( 1-|) 4b 2 x 2 

$ x) » 2 3> e s — 5 — 

o a 

o 


. (23) 


(24) 


This has the form of a harmonic oscillator potential 
with a frequency 


4 b 2 

, o° 2b (1-0 

or = 5— e 

m 


(25) 


from which y - b The w as a function of 

O 

volume can be found from equation (25) and the value 

of cu at £ = i. This can be obtained from entropy 
o 

measurements. In the limit of the Einstein model, 
the entropy becomes 


s = 3Nk (^^) - 3Nk ln (1 - e ' WkT) 

(26) 

-Nk (ln 8 - i) 

Given S and T, the above expression can be solved 
for cu. 


Figures 2 and 3 show how the vibrational com- 
ponent of heat capacity approaches ideal gas-like 
behavior in the limit of high temperatures or low 
densities. 


Electronic Contributions 


grand partition function, one obtains [9] 


N r e^ 2 de 


JN r 

V ® J 


0 e 


(e-M) kT 


and 


„ r e^de 

g j e (e-ft)/kT 


where 


2 l/2 Jrm z/2 


(27) 


(28) 


The chemical potential pi must be found for a 
given volume and temperature by a numerical iterative 
solution of equation (27) , then equation (28) may be 
integrated to find E. The pressure is obtained from 
2 

the identity, PV = — E, which holds for Fermi 

as well as Bose gases. Since, by definition, the 
chemical potential is the Gibbs energy per atom, 

G = N/^. Using the identity G = F + PV, the electronic 
contribution to the Helmholtz energy may be found. 


For temperatures that are small compared to 
the Fermi temperature, the integral in equation (27) 
may be approximated 


A 2 t 2 

12 e„ 


(29) 




and equation (28) becomes 


E = 


2 

5 


N € f + 


7I 2 

4 


N k 2 T 2 


where the Fermi energy e is 

r 





(30) 


(31) 


The first terms in equations (29) and (30) are the 
zero -degree degeneracy energy. This is already 
contained in the lattice terms and will therefore 
be suppressed. The electronic contribution to the 
Helmholtz energy is 


The electronic terms are obtained by treating 
the free electrons as an ideal Fermi gas. From the 


F = - 
e 


j^N k 2 T 2 
e 
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Figure 2. Behavior of heat capacity for solids and liquids in the high temperature regime. (Breaks 

in the curves represent phase changes.) 


The Fermi energy is computed from the observed 
electronic heat capacity coefficient which can be 
measured at low temperatures. The is obtained 

by differentiating equation (30) , giving 


COMPUTATIONAL RESULTS 
Specification of Constants 


c 


ve 



tt 2 Nk 2 T 
2e„ 


(32) 


At very low temperatures, the electronic terms 
dominate the heat capacity, which is experimentally 
observed to depend linearly on T. The value of 

obtained in this manner will be somewhat different 
F 

from the value obtained from equation (31) if the 
free electron mass is used for m. This difference 
arises from the interaction of the electrons with the 
lattice. The effective electron mass is defined as that 
mass which makes equation (31) agree with the Fermi 
energy obtained from the experimental heat capacity. 


The equation of state developed in the preceding 
section was applied to aluminum. The various 
constants required are obtained in the following 
manner. 

The zero-degree heat of vaporization L q is 

found from tables [ 10] to be 86 400 kcal/mole. The 
volume at 0°K is obtained by extrapolation of thermal 
expansion data and is taken to be 2.736 gm/cm 3 . 
From compressibility data, b is found to be 3.241. 

The Slater relation [11], 

2 1 

y = - — - — VP"/P' , where primes denote 

o 2 K K 
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Figure 3. Transition of heat capacity from liquid to gaseous behavior. (The curves indicate the behavior 
of the vibrational component of the liquid phase only, and ignore the liquid -gas coexistence region. ) 


differentiation with respect to V, gives y to within 
a few percent of the experimental value at normal 
density. Because of this agreement the Slater relation 
was used for the calculation for solid aluminum 
despite the objection that it assumes a volume 
independent Poisson ratio. In the development of the 
Slater relation, the result 


aluminum, the Einstein temperature, 9 = hco/k was 
found to be 269 , 38° K. For the liquid phase, the 
entropy in equation (26) is used to find the liquid 
Einstein temperature which is corrected to V q 

using equation (25) . For liquid aluminum, the 
Einstein temperature was found to be 117.935°K. 



V v > 

p k (V 0> 


is obtained. Given co , the Einstein frequency at 
o 

V , the 0 ) at any other V can be found. The 
o E 

quantity is found by comparing the entropy of an 

ideal solid in the Einstein approximation to the 

measured value at low temperature and using the 

Slater relation to correct to V . For solid 

o 


The electronic terms for aluminum are obtained 
by considering all three valence electrons as free and 
by taking the effective mass to be 1.6 times the normal 
electron mass. This brings the calculated electronic 
heat capacity, equation (32), into agreement with the 
measured value [ 12] . 


The constants Lm, B, E , £< in the liquid 
1 m 

potential are found in the following manner. At the 

melting volume ratio £ and temperature T , the 
m m 


1 
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energy is E (| )+ E (| , T ) + E (| , T ) 

K 111 v m m c m m 

which must equal to the energy of the zero -degree 

crystal plus the heat to melt, which is known 

experimentally. All quantities required to compute 

E (| , T ) and E (£ , T ) have been specified; 

v m e m m 

thus, is known. Similarly, 

P, (£ ) + P (£ , T ) + P (£ , T ) = P 
k s m v s m m' e s m m' atm 

therefore, P^Q^) is known. Using equations (15) 

and ( 16) for £ £*, simultaneous equations 

containing Iq and B are obtained. Their solution 
for aluminum yields Iq = 72 531 cal/mole, and 
B = 146 397 cal/mole. The quantity ^ is found by 
requiring the two expressions for P^, equation ( 16) , 

to be equal at £ = £ t . For aluminum, it was found 
that = 1.0265. Finally, the requirement that both 
expressions for E^, equation ( 18) , agree at 

£ = £i is used to obtain E = 2654 cal/mole. 

The Fusion Curve and Equation of 
State Surfaces 


Having specified the various constants for 
aluminum, computations were carried out for both 
the solid and liquid phase for various values of £ 
and T. Isotherms of a Gibbs energy versus pressure 
plot are shown in Figure 1. Since the phase that 
produces the lowest Gibbs energy at a given temper- 
ature and pressure is the stable phase, the inter- 
sections of the solid-liquid isotherms represent 
points on the fusion curve (Fig. 5) . Using the fusion 
curve, isotherms are constructed on a P-V surface 
(Fig. 6) , an E-V surface (Fig. 7) , and an S-V 
surface (Fig. 8) . Isoenergy lines on a P-V surface 
are shown in Figure 9. 

Comparison with Experimental Results 


To compare the equation of state developed 
in this paper with data obtained from shock com- 
pression experiments, the Hugoniot requirement, 

E - E = 4- P(V - V) , (33) 

o 2 oo 

is solved simultaneously with the equation of state. 

In experiments with porous samples, the volume 



Figure 4. Gibbs energy isotherms. (Circles 
represent the intersection of the liquid and solid 
isotherms and determine the melting pressure for 
that particular temperature.) 



Figure 5. Melting curve deduced from Gibbs 
energy plot. 

V qq represents the actual specific volume of the 
sample, which may be several times the volume of a 
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Figure 6. Pressure isotherms. (Breaks in the isotherms represent transition from solid 

to liquid phase. ) 


normal crystal. The states that can be reached 
by shocking a crystal with initial specific volume 

ratio V must lie along the intersection of 
oo 

equation (33) and the equation of state surface in 
P-V-E space (Fig. 9) . By varying the porosity of 
the sample, any state can, in principle, be reached 
by shock compression. 

In performing shock compression experiments, 
the observables are the shock velocity D and the 
material velocity u. These are related to the 
thermodynamical quantities E, P, and V through 
the Hugoniot relations 

V (D-u) = VD (conservation of mass) (34) 
oo 

PV = Du ( conservation of momentum) . (35) 
oo 

These relations are used to relate the observables 
u and D to the Hugoniot states found by solving 
equation (33) with the equation of state. Figure 10 


compares the predicted results for various porosities 
with measured results summarized in Table 1. Very 
good agreement is obtained considering that the 
equation of state constants are determined solely from 
ambient properties and contain no constants adjusted 
to fit the high-pressure data. The highest pressure 
datum point falls somewhat above the predicted curve 
for liquid aluminum, but is below the curve for solid 
aluminum. The fusion curve clearly indicates that 
melting should have occurred at this state, but the fact 
that the material was shock compressed may have 
prevented the phase transition to become complete in 
the short time involved, and the discrepancy could be 
that the aluminum is behaving as a super-heated 
solid. On the other hand, the inaccuracies in the 
model, particularly in the assumption of a Morse 
potential configured from the ambient compressibility, 
could easily account for such a discrepancy at 
pressures of 5 Mb. 

The reproduction of states on the Hugoniot for 
normal crystal density is not a particularly crucial 
test for an equation of state. The lattice terms are 
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Figure 10. D-u plots comparing computed results 
with experimental data points. (The short dashed 
region on the v oc / v 0 = i curve represents the 

transition region from solid to liquid phase. The 
long dashed curve represents a superheated solid. ) 

the major contributors at all but the highest com- 
pressions, and the fact that the restoring forces 
"stiffen up" at high compressions extends the validity 
of the assumption that the atoms vibrate as harmonic 
oscillators to high temperatures. Therefore, the 
Einstein or Debye Model together with the Mie- 
Gruneisen equation and a properly adjusted zero- 
degree isotherm will give a good representation of 
the normal Hugoniot to several megabars. Addition 
of electronic terms will extend such a model even 
further. However, such models fail for the lower 
compression, higher temperature states obtained by 
shock compressing porous samples. The equation of 
state developed in this work successfully predicts 
these states, as shown in Figure 10. 

It is useful for shock compression work to plot 
various thermodynamic coordinates against material 
velocity u as a means of expressing the thermo- 
dynamic state attained in a shock process. The 
material velocity is particularly useful because it can 


be easily related to the relative velocity of the 
impacting samples. Figure ii is a plot of pressure 
for different porosities versus u, and experimental 
points are shown for comparison. 

Figure 12 is a plot of temperature in the shocked 
region as a function of u. The various datum points 
are illustrated not as a comparison, since measure- 
ment of temperatures were not made, but as an 
indication of the range of temperatures accessible to 
this technique. The fact that higher temperatures 
are produced by increasing the porosity is dramat- 
ically illustrated by comparing the temperature for 
the 1.43, 2.08, and 2.98 points and 1.00 point at 
5.62 km/sec. These states were created by the same 
impact velocity. The fact that there is a discontinuity 



' w ( Iw/hc ) 

Figure 11. P-u plots comparing computed results 
with experimental data points. (No abrupt change is 
apparent that would indicate a phase transition.) 
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TABLE 1. SUMMARY OF HUGONIOT DATA FOR ALUMINUM 



GROUP I a 



D (km/sec) 

u^ (km/ sec) 

P(kb) 

v/v 

o 

7.531 

3.230 

335.8 

0.7874 

6.927 

2.319 

222.7 

0.8333 

6.500 

1.700 

153.5 

0.8696 


GROUP n b 



D (km/sec) 

u (km/ sec) 

P(kb) 

V /v 

O 

9. 13 

2.80 

693 

1.442 

10.39 

3.70 

1042 

1.553 

12.94 

5.62 

1971 

1.767 


GROUP 111° 


V 00/ V 300 

D (km/sec) 

P(Mb) 

V /v 

O 

d 

u 

1.00 

18.31 ± . 16 

4.93 

2. 185 

9. 93 

1.43 

11.74 ± . 10 

1.391 

1.498 

6.25 

2.08 

11.42 ± .09 

1.003 

1. 176 

6.74 

2.98 

10.75 t .08 

0.702 

1.015 

7. 18 


a 

These are representative points of those listed by J. M. Walsh, 

M. H. Rice, R. G. McQueen, and F. L. Yarger, Phys. Rev. ; 108 , 
196 ( 1957) The Material velocity u is taken to be i/2 the measured 
free surface velocity, u^. 

^L. V. Al’tshuler, S. B. Kormer, A. A. Bakanova, and R. F. Trunin; 
Soviet Physics — JETP, _ii, 573 ( 1962) 

°S. B. Kormer, A. I. Funtikov, V. D. Urlin, and A. N. Kiksnikova, 
Soviet Physics — JETP, _15, 477 ( 1962) 

^The quantity u was not stated in the paper. It was recovered by the 
relation u = [1 - (V/V ) (V 300 /V oo ) ] D. 



o 2 4 6 I 10. 12 14 

u(km/$tc ) 


Figure 12. T-u plots. (The data points are shown 
only to indicate range of temperatures reached by 
present experimental techniques. The break in the 
curves indicates that temperature measurement 
could serve as a method of detecting 
phase transitions. ) 

in the temperature curve at the melting region gives 
a possibility of experimentally detecting the melting 
point at very high pressures, although such temper- 
ature measurement would be extremely difficult to 
make in solids. 

Figure 13 is a plot of entropy as a function of u.. 
As before, the indicated data points are intended only 
to indicate range of experimental states, not measure- 
ments. For comparison, the temperature and phase 
of aluminum at ambient pressure for various values 
of entropy is indicated. Assuming the release is 
adiabatic and isentropic, an estimate of the release 
temperature can be made. Attempts have been made 
to measure the release temperature of shock 
compressed solids. Taylor [ 13] reports a favorable 



u ( km/ sec ) 

Figure 13. S-u plots. (The data points indicate 
measurement range rather than actual 
measurements. The dashed lines indicate entropy 
values corresponding to the stated temperatures at 
ambient pressure. If the release is assumed to be 
isentropic, these can be used to indicate release 
temperature. ) 

comparison between calculated and measured release 
temperatures for Cu up to the melting point. 

CONCLUSIONS 

A complete equation of state for metals that 
describes both the solid and liquid-dense vapor phase 
has been developed. All empirical constants required 
can be obtained from elementary thermodynamic data 
at ambient conditions . Since all thermodynamic 
functions are derived from the partition fimction, 
thermodynamic consistency is guaranteed. Phase 
transitions between the solid and liquid come about 
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naturally by using the Gibbs energy to determine the 
stable phase for a given state. Detailed computations 
were carried out for aluminum and the ability for the 
model to predict the behavior of aluminum for 
pressures of 5 Mb and temperatures of 20 000° K was 
demonstrated. Since only elementary thermodynamic 
data at ambient conditions are required to configure 
the model, extension to other metals is straight- 
forward by insertion of their appropriate constants. 


The model can be extended to higher temperatures 
with some additional computational effort to obtain 
the electron contributions for temperatures com- 
parable to the Fermi temperature. The model can 
also be extended to lower densities, provided the 
temperatures are low enough to completely neglect 
electronic contributions. 
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